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Abstract 

Charm dynamics is moving back into focus with the established discovery of oscillations in the neutral D meson 
system and the sign of direct CP asymmetry in D° — > tt + it~ /K + K~ . It opens the possibilities of finding CP 
violation beyond the reach of the Standard Model. In the recent past we have extensively studied charm dynamics 
within non-ad-hoc models with interesting flavour structures. We have shown that rare decays and CP violations 
are the best places to probe for New Dynamics in charm. We continue to study a different class of models, i.e., the 
Randall Sundrum model with a warped extras dimension to check for unusual effects in charm dynamics: namely in 
decays of final states with leptons and neutrinos and some asymmetries. These states should typically be accessible 
to experimental probes in the near future and, for certain, to any super flavour factory. 



Contents 



Contents 

1 Introduction 3 

2 Charm and the Standard Model 4 

2.1 D° - D°: Experiments and Theory 4 

2.2 Rare decays of charm 5 

2.3 CP violation in D° -> h+h~ 8 

3 A Randall-Sundrum model with custodial isospin 9 

3.1 Gauge structure and electroweak symmetry breaking 10 

3.2 Quarks and leptons 11 

3.3 Couplings of the fermion and gauge mass eigenstates 12 

4 The flavour structure and the parameter sets 14 

5 Rare Decays 15 

5.1 D° -» 77 16 

5.2 D° -> n+yr 16 

5.3 D -> X u vv 18 

5.4 D -> XJ+l- 19 

6 Asymmetries in D — > X u l+l~ 20 

7 Correlations between different AC = 1 observables 21 

8 Correlations between strange, charm and beauty 23 

9 Analysis of the parameter space 24 

10 Inference 25 

11 Summary 26 

12 Acknowledgement 27 
References 27 



2 



1 Introduction 



1 Introduction 

One puzzle in high energy physics is the hierarchy that exists between the scale of electroweak symmetry breaking 
(EWSB) and the Planck scale, Mpi, at which the effects of gravity must be taken into account. This Hierarchy 
Problem has inspired a great number of theories beyond the Standard Model (SM); one attractive solution being 
to assume the existence of a warped extra dimension. This class of models are known as Randall- Sundrum (RS) 
models [1]. They make use of the geometry of the extra dimension to provide a natural explanation for the hierarchy 
between the EWSB scale and Mp;. The RS geometry is built from a slice of AdS$ bounded by two branes, where 
SM particles are allowed to propagate in the bulk between the electroweak or 'infrared' (IR) brane and the Planck 
or 'ultraviolet' (UV) brane. In order to properly address the Hierarchy Problem, the Higgs is localized near or 
at the IR brane and the size of the extra dimension is adjusted to properly redshift the Planck scale towards the 
electroweak scale. In addition to providing an attractive solution to the Hierarchy Challenge, different incarnations 
of the RS model have been implemented to address the fermion mass hierarchy [2-5]. This is possible by adjusting 
the localizations of the fermion fields relative to the Higgs five dimensional background profile. The existence of 
a warped extra dimension and the ability of matter to propagate in the bulk leads to the existence of an infinite 
tower of Kaluza-Klein (KK) modes in the four dimensional effective theory. The existence of KK modes can lead to 
sizable contributions to flavour changing neutral currents (FCNC) well beyond limits given by the data. Variations 
of the RS scenario have been constructed in order to suppress these contributions and remain consistent within 
current experimental constraints [6,7]. 

The existence of electroweak precision data (EWPD) has placed very stringent bounds for the masses of KK 
modes. In particular, they constrain the masses of lowest KK mode excitations to Mkk > 12 TeV for a strictly 
IR localized Higgs field [8] and Mkk > 7 TeV for a Higgs field localized near the IR brane [9]. The LHC is on its 
way to probing large mass scales, but direct detection of new resonances will likely be possible close to the 1 TeV 
regime. This has motivated further research to extend the RS model in order to place KK ecitations within the 
reach of the LHC. In particular, work by the authors of [10] makes use of an extended gauge group in the slice of 
AdS§, protecting EW observables from very large contributions beyond experimental limits. They conclude that 
it is possible to fit EWPD with KK masses around 3 TeV. Furthermore, the authors in [11, 12] have argued that 
discovering these modes would be possible at the LHC with center of mass energies of 14 TeV and 100 fb _1 of 
integrated luminosities. Models that modify the AdS$ metric near the IR brane have also been shown to decrease 
the bounds on KK modes [9, 13-15]. Unlike the RS model, these extra dimensional structures are able to lower 
the KK mass scale to 2-3 TeV in large areas of their parameter space. Within this class of models, impact of KK 
modes involving third generation of quarks were shown to be within reach of the LHC with center of mass energy 
of 8 TeV and 10 fb _1 of integrated luminosities [16,17]. 

Dynamics of charm - including its mere existence - has historically been the proving grounds for many subtle 
aspects of the Standard Model. In spite of that, charm has been for long been treated as the step child in flavour 
dynamics for the lack of experimental statistics and theoretical handle over its dynamics 1 . The recent evidence for 
direct CP violation in AA C p by the LHCb [18], the CDF [19,20] and the Belle [21] collaborations has brought an 
end to those decades of neglect 2 . Although we do not address the implications of the class of models discussed in 
this work on the CP asymmetry in two-, three- and four- body hadronic decays of neutral charm mesons, and for 
very good reasons that we shall elucidate on later, we display other avenues through which warped geometry can 
leave its mark in charm dynamics. With the statistics available to the LHCb and in super flavour factories in the 
near future, these modes should be accessible through measurements of both decay rates and several asymmetries. 

In this work, we study the effects of KK modes to rare decays of neutral D mesons within and RS framework 
with custodial isospin protection. As explained above, this class of models leads to KK mode excitations which 
are currently being probed by the LHC or will be in the near future. The structure of this model is a lot more 
complex, with the existence of an extended fermion sector. Furthermore, the presence of KK gauge bosons leads 
to FCNCs already at the tree level which has interesting implications for CP violating observables as well as rare 
decays of K, B and D mesons. Signatures of this class of models along with the correlations generated by the same 
in K and B physics has been studied in great detail by [22-31]. They have extracted the flavour structure of these 
models as well as the fermion mass hierarchy which is consistent with current experimental bounds. Studies have 
also been done of the dependence of flavour dynamics on KK mass scales. In light of recent measurements by the 

1 "I know she invented fire - but what has she done recently?" 

2 Of course, we must not forget the "fans" of charm who supported it through out its not so glorious days and we pay them our due 
respect. CDF and Belle measurements are also sensitive to the individual asymmetries in D° — > K + K~ and D° —> ir + ir~ . 
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LHCb, CDF and Belle collaborations [18-21], we try to to compliment the work on K and B mesons combined 
with recent results on the B s — > // + /jT branching fraction [32], and study the rare decays of neutral D mesons. 

A new angle was introduced into our analysis by the very recent evidence of the time integrated branching 
fraction of B s fi + n~ at 3.5<r significance [32]: 

BR oxp (B s -> fi+fi-) = 3.2±\- 5 2 x 1(T 9 . (1) 

It is quite consistent with the time integrated SM expectation of [33] 

BR SM (B S -> M + lO = (3.54 ± 0.30) x 10~ 9 . (2) 

On the other hand this evidence leaves enough room for the intervention of new dynamics (ND) of several varieties, 
in this decay channel, it leaves us with the obvious doubt that SM might dominate over ND here and hence the 
latter might not be "found" here. We test our hypothesis of whether charm dynamics can allow significant ND 
intervention over SM effects even if the branching fraction of B s — > continues to look very akin to the SM 

expectation. While the SM produces the leading source of rare decays and asymmetries in B and K dynamics, 
the situation is very different for D transitions: SM gives tiny rates, and those are dominated by long distances 
dynamics. However, ND can enhance asymmetries of several classes significantly. Of course, one still has to deal 
with very small rates; however LHCb and super flavour factories will give us the needed huge data samples, while 
the SM might give little 'backgrounds'. 

The structure of this paper is as follows: In Section 2 we review the role of charm in the Standard Model and 
its current experimental status. We briefly discuss the different rare decay channels as well as the experimental 
evidence for CP violation in the D — > h + h~ channel at the Tevatron and the LHC. In Section 3, we introduce 
and review the main features of the RS model with custodial isospin protection. In Section 4, we discuses the 
flavour structure of the model along with the parameter space that we shall examine for AC = 1 processes. We 
calculate the new contributions to the effective Hamiltonian for the following AC = 1 observables in the decays: 
D° —> fJ. + n~, D — > X u vv and D — > X u l + l~ and discuss the results in Section 5. In Section 6 we look at the ND 
signature in several asymmetries in D — > X u l + l~ . We discuss the correlations between several flavour observables 
in Sections 7 and 8. In Sections 9, 10 and 11 we analyse and summarise our results. 



2 Charm and the Standard Model 

In this section we build up a case for studying the effects of ND in charm dynamics. We go over the SM estimates 
in multiple channels and the experimental observables that have been under scrutiny and some more which can be 
scrutinized in the near future to give us some insight into the basic dynamics that affect charm hadrons. 

2.1 D° — D°: Experiments and Theory 

That oscillations is a reality in neutral charm meson system is a well established result now [34-36]. The HFAG 
average stands at [37] 3 : 
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(3) 



If it is assumed that there is no direct CP violation in charm dynamics, values of \q/p\ and (f>o come closer to their 
benchmark values for no CP violation: 



= 1.04±83S. ^=(-2.02t^)°. 



(4) 



Recently, the LHCb Collaboration has presented measurements of oscillation parameters and ruled out the no 
oscillation hypothesis at 9.1a [38]. This is the first exclusion of the no oscillation hypothesis made from a single 



Up to date results can be found in the Heavy Flavour Averaging Group's website. 



4 



2 Charm and the Standard Model 



measurement. The measurement of the time dependent ratio between the doubly Cabbibo suppressed D° — > K + ir 
to the Cabbibo allowed D° —> K~n + yields: 

x' 2 = (-0.9 ± 1.3) x 1CT 4 , y' = (7.2 ±2.4) x 10~ 3 , R D = (3.5 ± 0.15) x 10~ 3 . (5) 

However, the relative size of xd and yo is not clear yet. Before these experimental results in 2007 [34-37], most 
authors had argued that the Standard Model (SM) predicts Xd, yrj < 3 x 10~ 4 — yet not all. In 1998, xd, yo 
< 10 -2 was listed, admittedly as a conservative SM bound [39], together with a question: How can one rule out 
that the SM can not produce 10~ 6 < r D < 10~ 4 (corresponding to x D , y D ~ 10~ 3 - 10" 2 ). In 2000 and 2003, 
an SM prediction obtained from an operator product expansion (OPE) yielded xd, yD ~ O(10~ 3 ) [40,41]; later 
another OPE analysis was done with similar results [42,43]. Alternatively in 2001 and 2004, an SM prediction on 
D° — D° oscillations was based on SU (3) breaking mostly in the phase space for yo and then from a dispersion 
relation for x D [44,45]. 

The question of whether xn and yr) can be accommodated within the SM is still open. Furthermore the 
situation is more complex: many theorists think that yo depends almost all on long-distance dynamics and is 
hardly sensitive to impact from ND - in contrast to xd that can be produced by ND (unless sources of ND suffer 
cancellations). 

While hints of CP violation have been seen in recent experiments, it is still not clear whether these can be 
interpreted as a distinct signature of new dynamics (ND). However, there have been suggestions in the past that 
ND can possibly make large impact in charm dynamics [46-48] . The effect of models with a warped extra dimension 
in D° — D° oscillations has been considered in [28] and it was shown that certain parts of the parameter space can 
be constrained by the experimental limits in these observables available at that time. 

2.2 Rare decays of charm 

Rare charm decays within the SM have been studied extensively in the literature. The verdict is universal: SM 
leaves very small effects in these rare decays. Table 1 summarizes the SM effects in the rare decay channels that 
we will look into in this work. The results can be summarized as follows: 



OBSERVABLE 


SM SD 


SM LD 


EXPERIMENT 


BR(£>° - 


-> 77) 


(3.6-8.1) x 10~ 12 t 


(1 - 3) x IO" 8 [49,50] 


< 2.4 x 10~ 6 [51] 


BR(£>° - 




6 x 10~ 19t 


(2.7-8) x 10~ 13 


< 1.3 x 10~ 8 


BR(D - 


> X u vv) 


10~ 15 - io- 16t 






BR(D ± 




3.7 x 10~ 9t 


- o(io- 6 ) 


- o(io- 5 ) 


A c 

^FB 




- 2 x 10~ 6t 










- 3 x 10~ 4 f 






A CP 
^FB 




- 3 x 10~ 5t 







Table 1: SM short distance (SD) and long distance (LD) contributions to D° — > 77, D° —> , D — s- X u vv and 

D ± — > X u l + l~ . Detailed calculations of the numbers marked with a (t) can be found in the references [52-54]. 



D° ->■ 77 

The short distance contribution to this channel is about three orders of magnitude smaller than the long 
distance contribution. The long distance contribution comes primarily from a vector meson intermediate 
state and has been calculated in [49]. Alternatively, the long distance contribution has also been estimated 
using the hybrid chiral perturbation theory in [50]. New dynamics can possibly contribute to this channel 
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through short distance operators and hence has to enhance the short distance SM contribution by at least 
three orders of magnitude in order to overcome the SM LD contribution and make its presence felt. In our 
previous work on LHT-like models [52], we showed that new dynamics of this class can enhance the short 
distance contribution but not enough to overcome the long distance one. 

Current experimental limits on this channel are quite far from the the SM estimates. Recent results from 
BABAR cite [51]: 

BR oxp (L>° -> 77) < 2.4 x 1(T 6 at 90% C.L. (6) 
BESIII has updated their results recently [55,56]. However, their limit falls short of the limit set by BABAR 

BR B esiii(£>° 77) < 4.6 x 1(T 6 at 90% C.L. (7) 

The preliminary BESIII result is with 2.9fb~ 1 data at the ^"(3770) peak. They have a projected reach 
of 5 x 10~ 8 with 10fb _1 at the i/>"(3770) [55] peak within the next few years. A good measurement of 
the branching fraction of this channel will give us a very good estimate of the SM LD contribution to the 
branching fraction of D — > as explained below. 

D° -»■ 

The short distance contribution from SM to this channel is extremely tiny and is some six orders of magnitude 
smaller than the SM long distance contribution. The latter is driven by the two photon unitary contribution 
[49] and hence depends directly on the size of BR(Z?° — > 77) through the relation 

BR^(L>° fi+fi-) ~ 2.7 x 10~ 5 BR(£>° ->■ 77). (8) 

New dynamics that can enhance BR(£>° — > 77) above its SM value will indirectly enhance the SM LD 
contribution to BR(D° — > /i + /i~). It is also believed that any new dynamics that contributes to D° — D° 
oscillations will also contribute to this decay channel significantly [57]. This is in general true. However, it 
should be kept in mind that such new dynamics can enhance the short distance SM contributions by many 
orders of magnitude but might fail to overcome the long distance contribution as we showed in [52]. The SM 
SD contribution to the branching fraction of D° — s- n + ijT is given by [58]: 



BR(.D - fyr) = (^^) 2 

x E KjVuj^Yoix^ + ^ixj)) 2 , (9) 

j—d : s : b 

where — 11 2 /m^. The definitions of Y (xj) and Y\ (xj) can be found in our previous work [52]. ND from 
the warped extra dimension can enter through tree level contributions to Y(x) as defined later in Section 5. 

The LHCb now reports an upper bound of [59]: 

BR cxp (L>° ->■ fi+fj,-) < 1.3 x 10~ 8 at 95% C.L.. (10) 

with 0.9fb _1 of data. They can be expected to enhance this measurement by two or three orders of magnitude 
in the future. This measurement is out of reach for BESIII and can be within reach of super flavour factories. 

D — > X u vD 

The SM contribution to this channel is extremely small. A naive estimate using the Inami-Lim functions [60] 
leads to the numbers [49]: 

BR SD ( J D+ -> X u vv) ~ 1.2 x 10~ 15 , BR SD (L>°^X u ^)~5x 10" 16 . (11) 

A complete analysis of SM contributions to this channel requires a full two-loop renormalization group analysis 
[61] as it is done for the light quark loops in the analogous K decay channel. This is necessitated by the fact 
that the quarks running in the loops of the penguin diagrams are light states and hence these operators cannot 



6 



2 Charm and the Standard Model 



be treated as purely short distance ones. Alternatively, these contributions can be treated as vector mesons 
in the intermediate state and such long distance contributions can be estimated [49]. The enhancement of 
such contributions over the SM SD contribution can be about an order of magnitude. 

The impact of ND on this channel can be significant. In our analysis of this channel within the LHT-like 
framework [54] we showed that enhancements of O(10 3 ) — O(10 4 ) are possible. Considering the experimental 
sensitivity attainable in the immediate future, attaining theoretical precision in the calculation of the SM 
contribution to this decay channel is not a very fruitful enterprise. In spite of this, these channels are very 
demonstrative of the impact of ND over short distance SM dynamics and in the future super flavour factories, 
or in any future super charm factories, these decay channels with the requisite ND enhancement might open 
up for experimental probes. 

D ± ->■ XJ+l- 

In both the inclusive and exclusive modes, the branching fraction is dominated by long distance effects driven 
by the p, u and <f> mesons. In size they are larger by two or three orders of magnitude than the quark level 
contribution to the inclusive mode, depending on the exclusive mode being looked into. The short distance 
contribution is primarily driven by the photon penguin and new dynamics of the LHT-like class can hardly 
make a dent there [53]. However, one can always look at the dilepton invariant mass distribution of the 
differential decay rate and distinguish between the short and the long distance contributions. Any deviation 
from the shape or size of this distribution as predicted by the standard model will be an unambiguous signal 
of new dynamics. The SM SD differential branching fraction for the inclusive process D — > X u l + l~ is given 
by [53]: 



^BR^ (D -> XJ+l- 

1 G\o?m\ 3 
T D 768vr 5 1 



(|C 9 ( M )| 2 + \CM\ 2 ) (1 + 2S) + 12 Re(0 7 (/z)C 9 *(/i)) +4 ( 1 + - J \C 7 (ja)\ 



(12) 



with 



(Pi+ + Pi- f 

S = n 



We set fj, = m c = 1.2 GeV. Integrating over s gives us the total decay rate. One has to be careful about 
not picking up the infrared divergence in the differential decay rate. We made an infrared cut on s at about 
an invariant dilepton momentum of 20 MeV. The definitions of the Wilson operators and the form of the 
coefficients can be found in [53] . Models with a warped extra dimension can leave their impact in Cg and C\o 
trough tree level contributions to Y(x) and Z(x) defined below in Section 5. 

Asymmetries in D^ 1 — > X u l+l~ 

Although the branching ratios are dominated by long distance physics, the asymmetries are sensitive to mostly 
short distance physics. We clearly showed in [53] that if one kinematically cuts out the vector resonances 
in the invariant dileptonic mass distribution, one loses effects of 0(10%) only in the asymmetries. Hence 
asymmetries are good observables for the discovery of ND as the latter is expected to bring enhancements 
to the short distance dynamics. There are three asymmetries that can be probed in this channel which come 
solely from AC = 1 currents [53]. 

The forward-backward asymmetry Ap B is tiny in the SM, O(10 -6 ), and the impact of new dynamics here 
can be quite large. While this held true in the LHT-like model we looked at in our previous work [53], it 
was also true that the absolute size of the asymmetry remained at 0(1%). The normalized forward-backward 
asymmetry is defined from the double differential decay rate as 

/U^r(D±^XJ+|-)]sgn(*)ds 
^pr(s) = 1 5 ■ (13) 
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After performing the integral over the angular distribution we obtain 

-3 mC*MC 9 (ri)s + 23*(C 1 * ( M )C 7 (m))] 



^fb(s) 



(1 + 2s) (iCg^)! 2 + |C 10 (m)| 2 ) + 4 \C 7 (n)\ 2 (1 + f ) + 125i (C 7 ( M )C|(m)) ' 

(14) 



Enhancements to both Cg and Cio can manifest themselves as sizable A FB . 

The CP asymmetry A^ p is of O(10~ 4 ) in the SM and new dynamics like the LHT-like models can hardly leave 
much dent there. This can be well understood in terms of the absence of tree level contribution to the CP 
asymmetry in these class of models. However, this is not true for the models with a warped extra dimension 
and hence we check for enhancement to the CP asymmetry. The CP asymmetry parameter Aqp(s) is defined 
as 

^T(D+ -> xj+i-) - £t(d~ -> xj+i-) 
CPW ±T(D+ -> xj+i-) + ±T(D- -> X a l+l~) ' 

Integrating over the invariant dileptonic mass we get the total CP asymmetry A^ p 

The CP asymmetry in the forward-backward asymmetry A FB can show very large contributions from new 
dynamics like the LHT. The SM contribution to this asymmetry is of O(10~ 5 ). New dynamics can enhance this 
by an order of 0(10% — 100%) [53]. Since this asymmetry is sensitive to phases in the Wilson coefficients, it is 
open to enhancements by ND, and there is a good possibility that models with a warped extra dimension can 
enhance this asymmetry to quite large values. The normalized difference in the forward-backward asymmetry 
in D — > X u l + l~ and D -> X a l + l~ is defined as [62] 

^fb(s) = + ^tli^i (16) 

^fb(^) — A FB (s) 

In the limit of CP symmetry, A FB (s) and A FB (s) have to be exactly equal in magnitude but with an opposite 
sign [63,64]. As the forward-backward asymmetry is defined in terms of the positive anti-lepton, A FB (s) and 
A FB (s) have opposite signs. ^4pg(s) is sensitive to the phase in C7, Cg and C w . The SM offers phases only 
in C7 and Cg in D — > X u l + l~ and none in Cio- Hence the integrated asymmetry turns out to be very small. 

i4gg(S)dS = AgI~3x HT 5 . (17) 

To make the study of these asymmetries "clean" both theoretically and experimentally, cuts in the dilepton 
mass distribution can be made around the p, lo and (f> resonances. Our previous study showed that making 
such cuts left our results unaltered proving that long distance dynamics, especially the resonances, have very 
little to do in these asymmetries. It should also be emphasized that the procedure of kinematically "cutting" 
out the resonances can be done both in theory and experiments leaving equivalent impacts on the observables. 

2.3 CP violation in D° -)■ h+h~ 

There are two interesting results for CP violation in D° — > Tr + rr~ / K + K ~, both from hadronic experiments. The 
CDF results can be summarized as [19]: 

(A^ F (D° -> tt+O) = (+0.22±0.24 stQ i±0.11 syst )%, 
{A^ F {D° -> K+K-)) = (-0.24 ± 0.22 stat ± Q.lQ syst )% . (18) 

The result from LHCb [18] is cast in a different form to cancel inherent production asymmetries: 

AA^> Ch = Acp{D° -> K+K-) - A CP {D° -> tt+tO = [-0.82 ± 0.21(stat.) ± O.U(sys.)} %. (19) 

CDF also quotes AA CP [20]: 

AA^ F = [-0.62 ± 0.21(stat.) ± 0.10(sys.)] %. (20) 
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Of course, we should not overlook the recent Belle result using 976fb 1 of data [21]: 




(21) 
(22) 



The CDF and Belle results are sensitive to both direct and indirect CP violation. However, the difference in the 
CP asymmetry parameters, AAcp measured by all three experiments is sensitive mostly to direct CP violation. 
This holds true as it is well agreed upon that indirect CP violation has to be universal and hence mostly final state 
independent, which leads to a cancellation when differences in CP asymmetry in two modes are taken. It should 
be kept in mind though that AAcp does contain a part of the indirect CP violation through the interference of 
mixing with direct CP violation. 

A second subtlety that should be kept in mind is that it is generally agreed upon that the sign of the direct 
CP asymmetry in the K + K~ decay mode is opposite to that in the ir + ir~ decay mode. However, the sizes of the 
asymmetry need not be related and are equal only in the limit of exact SU (3) p symmetry. 4 This means that the 
asymmetries really add up in the difference. 

For indirect CP violation, significant deviation of the observed values of \q/p\ from 1 and 4>d from or 7r opens 
up possibilities of large CP violation. On the other hand, tiny values for xp> and yp, effectively suppress these 
effects. Hence, the fact that indirect CP violation is expected to be small in charm mesons is a consequence of the 
size of the oscillation parameters and is not related to the phase contributions to CP violation per se. 

Within the SM, it is believed that direct CP violation in these modes should be of O(10 -4 ) and indirect 
CP violation should be of O(10~ 5 ) [66-69]. Direct CP violation of O(10~ 3 ), however, cannot be claimed as an 
unambiguous signal of new dynamics. For that claim to hold, direct CP violation of O(10~ 2 ) has to be observed. In 
quite a few scenarios of new dynamics indirect CP violation can be enhanced without enhancing direct CP violation 
much above what the SM has to offer. Motivated by the initial results of CDF, we had shown that scenarios like 
the LHT cannot enhance direct CP violation while it can saturate the limits of indirect CP violation [66]. There 
are some others in which both can be enhanced. A lot of work has been come forward in the literature since the 
announcement of the LHCb results showing how direct CPV can be significantly enhanced by different scenarios 
of ND [6,70-81]. However, it can be equally well argued that the SM can accommodate the present experimental 
evidence [70,82-85]. Considering the uncertainties shrouding the SM estimate of CPV in this channel and moreover 
the conceptual shortcomings of the theoretical technologies used to make these estimate, it is quite impossible to 
make a simple statement of whether the measured value of Aj4cp is within the SM reach or is the signature of 



The work which is most closely related to the kind of models we will study in this work is reference [77], 
where an estimate is made on how the localisation of the Higgs in the bulk impacts the direct CP asymmetry 
in D° —> n + ir~ / K + K~ . In the light of the conceptual uncertainties plaguing the precise determination of both 
the SM and the ND contribution to the direct CP asymmetry in hadronic decays of charm, we shall refrain from 
making an estimate of the same. Suffice it to mention here that it can be inferred that models of this kind can well 
produce the CP asymmetry under question. Whether that is what is actually happening can only be answered if 
and when more rigorous theoretical tools have been developed and/or experiments definitively make a case for SM 



3 A Randall-Sundrum model with custodial isospin 

In this section we review the RS model with custodial isospin protection [10]. We examine how electro weak 
symmetry breaking (EWSB) proceeds and its implication for the spectrum, in particular how it leads to FCNC 
at tree level. Additionally, we review the larger fermion structure of this model and how it leads to fermion mass 
mixing beyond the CKM structure. 

4 More specifically, the asymmetries are equal in magnitude and have opposite sign in the limit of exact U-spin symmetry. The 
implications of U-spin was first pointed out in the context of DCS decays Gronau and Rosner in [65] and later in the context of CP 
violation in SCS decays by Grossman et. al. in [69]. We are assuming here that isospin symmetry holds good. 



ND [70,71,82,86]. 



or ND. 
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3.1 Gauge structure and electroweak symmetry breaking 

The RS model with custodial isospin consists of an extended gauge symmetry in the bulk given by: 

SU{3) C x SU(2) L x SU(2) R x U(l) x x Plr- (23) 

The model is constructed in a slice of AdS§ described by the following metric: 

ds 2 = e-^rj^dx^dx" ~ dy 2 . (24) 

The UV brane is located at L = and the the IR brane at a distance y = L along the extra dimension. In the 
bulk, the Lagrangian for the gauge sector is given by 

L gauge = s/G (^TtW M nW mn - \ttW M nW mn - \b mn B mn - \ttG mn G mn ^ , (25) 

where Wmn is the field strength for the SU(2)l gauge group, Wmn for SU{2)r, and Bmn and Gmn for U(l)x 
and QCD gauge groups respectively. By assigning the appropriate boundary conditions, (UV,IR), one is able to 
recover the electroweak (EW) gauge symmetry in the UV brane: 

W^ 2 < 3 (++), £„(++), 
W2' 2 (-+)> W*(++), (26) 

where (+)/(—) denote Neumann/Dirichlet boundary conditions. Furthermore, breaking of the bulk symmetry in 
the UV brane down to SU(2)l x U(1)y leads to mixing between £? M and W»: 

z l = Wl cos 4>- B^ sin <j>, 

B^ = sin </) + Bp cos <j>, (27) 



where 



COS0 = 

VW+M 



95 

^2 : 



95 (28) 



This symmetry breaking pattern guarantees massless zero modes for the fields with (+,+) boundary conditions. 

Electroweak symmetry breaking is achieved as in the RS model [1], introducing a scalar Higgs field, which is 
now localized near the IR brane. This field which transforms as a bidoublet of SU(2)r x SU(2)l is given by 



TT^ h°-i 

V2 2 

h °+iTT 2 TT- 

2 72 



(29) 



The neutral component of the Higgs field acquires a vev which breaks the EW bulk symmetry, SU(2)r x SU(2)l, 
down to its diagonal combination SU(2), leading to an unbroken custodial symmetry. As can be seen from the 
pattern of EWSB in the UV as well as the IR branes, the discrete symmetry, Plr, remains unbroken, essentially 
protecting non-oblique parameters from large contributions beyond experimental limits. 

After EWSB, the photon and gluon zero modes remain massless. Zero modes of fields coupling to the Higgs, 
W^' 2,3 , acquire a mass and the following mixing pattern arises between the neutral components of the gauge fields: 

z n = w l cos ^ _ B » sin 

A,, = Wl sin %j} + B^ cos ip. (30) 
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The mixing angle tp can be expressed in terms of the angle 4> using the following relations: 

cos ip 



y/l + sin 2 (f> 
VI + sm 2 <f> 

Within this framework the physical Z mass is generated through mixing between the zero and first 5 KK modes of 
Zfj, with Z' ; similarly for the charged gauge bosons. 

As in the RS model, profiles for the gauge boson KK modes are given by 

where J\ and Y\ are Bessel functions of the first and second kind, m, n is the n th order KK mode mass, and N n is 
a normalization factor. 

3.2 Quarks and leptons 

Fermions in this model have to be embedded in representations of the bulk symmetry since they propagate in the 
bulk. Additionally, contributions to EWPD beyond experimental limits need to be suppressed, such as Zbb. This 
can be done by embedding fermions within representations of the group 0(4), which is isomorphic to SU(2)l x 
SU(2)n x Plr [87]. As a consequence, this embedding will yield a larger fermion structure. 

In particular, left handed up- and down type quarks are embedded within bidoublet representations of S(2)l x 
SU(2)r, with the right handed up type quarks transforming as singlets under 0(4). The right handed down type 
quarks transform as a (3,1) ® (1,3). As with gauge bosons, fields with Neumann boundary conditions at both 
branes, (++), give rise to massless zero modes, leading to the SM model spectrum. The electric charges are given 

by 

Q = Tl + T% + Q x , (33) 

where Tj* R denote the third component of SU{2)l,r and Qx the U{\)x charge. Through this framework, all 
fermions receive the same U(l)x charge and after electroweak symmetry breaking, mixing occurs between fermions 
with the same electric charge. A detailed study of the mixing between zero modes and the lowest KK modes is 
carried out in [22]. Here we just summarize the most important results that lead to the quark mass matrices. The 
structure of the model consists of three mass matrices for the states with charges 5/3, 2/3, and —1/3. For states 
with left handed chirality, these can be grouped into the following vectors: 



*l(5/3) t = (xT(-+),i>'L(+~),K(+-)) Qx=2/3 , 



* L (2/3f = (# (0) (++),# {++\U'l{+-),U'l i {+-),xi{-+),<{— )) 
*l(-1/3) T - (qt ii0 \++),qt i {++),D'i(+-),D i L (--)) 



2x=2/3 



Qx=2/3 

(34) 



Similarly, right handed fermions are obtained by replacing L — > R in the above expressions. 

The structure of the Yukawa interactions have to preserve the 0(4) symmetry. This leads to the following four 
dimensional Yukawa matrices: 

(^).,-VSiSjf*»5^(»)^(»W»). 

5 We do not consider the higher KK modes as they have progressively smaller impacts. 
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where the functions Jl,r are the fermion 5D profiles given by [2,4] 
for the zeroeth mode and 



Hr (y) = 




with a = | ± c + 1/2| for the n th KK mode. X u ' d denote the fundamental 5D Yukawa couplings and c the bulk 
mass of the 5D fermion field. The Higgs field background profile is given by the function h(y). Using the weak 
eigenstates in Equation (34) and the effective 4D Yukawa couplings in Equation (35) one can write the following 
Lagrangian which mixes zero and first KK modes with the same electric charge [22] 

C = -*i(5/3)^t(5/3)*ji(5/3) + h.c. 

- ^ (2/3)7W (2/3) 4^(2/3) + h.c. 

- # L (-l/3)M(-l/3)* fl (-l/3)+ h.c. (38) 

The lepton sector can be accommodated in a similar fashion, but with the appropriate choice ofU(l)x charge, Qx —> 
0. Within this model, left handed charged leptons and neutrinos are grouped within a bidoublet of SU (2) l x SU (2)r, 
and the right handed charged leptons transform as a (3, 1) © (1, 3) of the EW symmetry. In addition, this model 
contains a right handed neutrino field transforming as a singlet of the 0(4) x U(V)x bulk symmetry. 



3.3 Couplings of the fermion and gauge mass eigenstates 

To calculate the tree level contributions to FCNC's we first derive the couplings of fermions to the electroweak 
neutral gauge bosons. Because of the mixing that arises between the zero and first KK modes of with Z' , we 
must diagonalize the mass terms in order to find the appropriate mass eigenstates. In this section we will work 
in the EW basis for fermions. The rotation to mass eigenstates will be shown in the following section, where the 
CKM structure of the model will be discussed in greater depth. In order to derive the couplings between neutral 
gauge bosons and up type fermions we follow the analysis in [87]. The left handed up type quarks transform as a 
bidoublet of the full 0(4) x U(i)x symmetry and couple to the gauge bosons through the following Lagrangian: 

C D Cl Tr[Q L ^Q L %] + c 2 Tr[Q Ll ^Q L ] + (c 3 Tr[Q L ^iD^U]Tr[U^ Q L ] + h.c) , (39) 

where Q L G (2,2) 2/3 of 0(4) x U(l) x , V„ = {iDJJ)U\ and % = {iD^lTjW. We make use of the following 
covariant derivative 

D^U = d^U + i 9 -^<J a W;U + i 9 -^a a W;U - B„Ua 3 , (40) 

where U is a non-linear sigma field which contains the pseudo Nambu-Goldstone bosons resulting from the breaking 
of 0(4) into 0(3). The Plr symmetry of the model leads to the conditions 95 = ^5 and c\ = C2; guaranteeing 
the vanishing of a potentially large contribution to ZbjJjL coupling which would violate experimental limits. We 
have naively used effective four dimensional fields together with the fundamental 5D gauge coupling. To obtain the 
effective 4D coupling, integration over the fermion and gauge boson profiles must be carried out with U = 1-2x2 
after the 0(4) symmetry is broken. Contributions to the right handed up type quarks vanish, since these transform 
as singlets under the 0(4) symmetry, coupling in the usual sense to the U(l)x gauge boson. 
The couplings of left handed up type quarks to the Z°, Z 1 and Z' are given by 

9l,z> = — A \-\ - \^ 2 <P] F' L (cL,i,CL,i), (41) 
cos0 \ 2 6 / 
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and those of right handed up type quarks by 

2 54 



9r,z° 



3 cos ^ 



sin 2 ipF R (c Rti ,c Rti ), 



2 54 .9 /, r i / \ 

5fl,z! = -g T sm vF R {c Ri ,c Rii ), 

3 cos 

5h,Z' = -|-^r sin2 <t>F' R (cR,i,CR,i)i ( 42 ) 

o COS y) 

where 54 denotes the SU{2)l coupling strength as in the SM model, and F' LR ' (cl,ra, cl.ua) are diagonal matrices 
that quantify the overlap between right and left handed zero-mode fermions with the Z' 0,1 gauge bosons. These 
matrices give rise to FCNC's at the tree level given that they are not proportional to the identity. The coupling to 
the photon is given by 

2 

9l,r = -gmmip. (43) 

The overlap is identical to that of the Z 1 gauge boson since at this stage they both have the same 5D profile. The 
effects of EWSB arise at the mixing level between the three neutral gauge bosons. The mixing between the neutral 
gauge bosons can be parametrized by a unitary matrix GZ used to diagonalize the masses of the neutral gauge 
bosons [22]. Therefore, the mass eigenstates can be ordered in the following way: 

(Z°,Z\Z') T = GZ(Z 3 ,Z 2 ,Z 1 ) T , (44) 

where Z\ denotes the lightest eigenstate and the couplings in Equations (41) and (42) transform such that gL,R,{i) ~ > 
9L,R{i)GZij. 

The analysis of the interactions between neutral gauge bosons and leptons is carried out in the same way. 
The interactions of the left handed leptons are obtained from Equation (39), which leads to the following effective 
couplings: 

54 

9l,z« = 



Qsin 2 ^J Fl(c L ,i,c LA ), 
g L , z >= -^(\sin 2 ^Ftic^c^), 



94 ( 1 , 1 _. 2 



9l,z> = -^-(-- + -sin 2 ^)F' L (c Lti ,c Lti ). (45) 

COS (f> \ I I J 

Right handed leptons belong to the (1, 3) triplet representation of SU(2)l X SU(2) R . An effective Lagrangian anal- 
ysis yields two additional contributions to their couplings to neutral gauge bosons [87] . These can be parametrized 
by the following effective Lagrangian 

C D c 4 Tr[E R YE R V^} + c b Tr[E Rl ^E R ] . (46) 

The effective couplings are given by: 

9r,z° = -^-7 sin 2 tl>F%(c Ri ,c Rti ), 
cosV' 

54 . 2 / Ell I \ 

9r.z^ = rsm ipF R {cR^,c RA ), 

cost/ 1 

9r,z< = -9acos())F r {c R a,c R a), (47) 

while those to the photon KK mode are given by 9l,r = —g4,shiip. 

In Section 5 we show the new tree-level contributions to the rare decays of D mesons due to the couplings 
introduced above using a KK mass scale of 2.45 TeV. Before proceeding, in Section 4 we discuss the flavour 
structure of the model and the parameter sets used in the analysis of D meson rare decays. 
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4 The flavour structure and the parameter sets 

The flavour structure of this model was studied in depth in [22,24]. Alternative approaches can be found in [28,88]. 
Before looking at the structure of the CKM matrix, it is worth emphasizing that the mixing of zero modes with 
heavier KK modes, Equation (38), leads also to a source of FCNC. However as it was pointed out in [24] as well 
as [6,88], this effects are negligible and thus we only consider the k = I = components of the 4D Yukawa matrices, 
Equation (35), when rotating to the mass eigenstate basis. The 4D Yukawas can then be effectively written as 

kL 

V/ = Kf-^Iliv = L^ Q )f Q R {y = A< d ), (48) 

where X"j d denote the 5D Yukawa couplings. The CKM matrix is obtained as in the SM, that is 

Vckm = U* L D L , (49) 

where Ul and Dl are unitary matrices which rotate flavour eigenstates into mass eigenstates for left handed up- 
and down type quarks. Their complete parametrization can be found in [22,24]. Since we are considering neutral 
current exchanges at tree level, we found it useful to define effective unitary matrices for the up type quarks, f/r,, e // 
and Uit eff coupling quark mass eigenstates to the three physical massive neutral gauge bosons as well as to the 
first KK mode of the photon. These can be obtained using Equations (41), (42), and (43): 

UL,eff(n) = g L ,z° ■ GZ ln + ulg L}Z iU L ■ GZ 2n + W L g L! z'U L ■ GZ 3n , 

Un,eff(n) = g R , z ° ■ GZi n + U R g R . z iU R ■ GZ 2n + U R g RjZ >U R ■ GZ 3n , (50) 
for the three physical neutral gauge bosons denoted by n — 1, 2, 3 and 

U L ,eff{i) = Ulg L U L , 

U R , ef fh') = U R g R U R , (51) 

for the first KK mode of the photon. In what follows we describe how the parameters encoded in the matrices gL,R 
are constrained, in particular, their dependence on the bulk mass parameter. 




Figure 1: Scatter plot of the Yukawa couplings of the quarks. The figure on the left is for the up type quarks and 
that on the right is for the down type quarks. 

Being motivated by "naturalness" and an attempt to generate a natural hierarchy from the Plank scale to the 
electroweak scale, the models with a warped extra dimension are at their best when the absence of fine tuning 
is prevalent in all parameters. However, severe constraints from flavour dynamics, especially from the measured 
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Figure 2: Scatter plot of the bulk mass parameter of the quarks. The figure on the left is for the left handed quark 
doublets and that on the right is for the right handed up type quarks. 

value of ex force some degree of fine tuning in these models. Even if these models have a custodial symmetry that 
protects them from electroweak precision constraints, further constraints on the KK mass scale are generated due 
to constraints from €k ■ However all this can be alleviated by a judicious choice of the parameter space that we 
deem as allowed. A "judicious" choice always means the reintroduction of tuning of the parameter sets. However, 
what is important here is not the presence or absence of tuning, but the degree to which this is necessary. The 
primary aim of this is to keep the KK mass scale as low as possible to make the model testable in the near future, 
while not allowing the violation of flavour constraints. 

In the light of this argument we shall display the parameter set that has been used in the analysis in a manner 
that lays it open to naturalness arguments. Amongst, the parameters relevant to this analysis only a subset yield to 
this argument on naturalness, namely, the Yukawa couplings of the quark sector. We do not discuss the naturalness 
in the Yukawa couplings of the lepton sector as it is not relevant to this analysis. 6 

The Yukawa coupling of the up type quarks seem to be quite "anarchic" and random as can be seen from 
Figure 1, being evenly distributed over the entire parameter space. On the other hand, the down type Yukawa 
couplings seem to have a tendency to cluster towards the higher values of the couplings for the first and the second 
family. It must be kept in mind that the allowed parameter set is shaped by both constraints from electroweak 
precision tests and flavour constraints, some of which are stronger on the first two families than on the third, 
especially for the down type quarks. 

In addition to the Yukawa couplings discussed above, it is also informative to take a look at the quark bulk 
mass parameters, as an interesting pattern is seen when those of the first two families are compared to the third 
family, Figure 2. For the quark doublet mass parameters CQ i: the allowed parameter space is a cone starting from 
high values of cq 1 . For the right handed bulk mass parameters c UjCi t, the cone is less tight and starts from lower 
values of c MiC; t. The very large positive values of c t should be noted, these are required to generate a the right value 
of the top mass, and thus large overlap between tL, tn with the Higgs profile is necessary. 

5 Rare Decays 

In this section we discuss the effects of the class of models introduced above on certain rare decay channels of 
charm mesons. As should be evident from the previous sections, FCNC from this model of ND comes from tree 
level exchanges of neutral gauge bosons. The primary contribution on all the decay channels that we have studied, 
and elucidated on below, comes from the mixing of the new gauge boson states, into Z\, the SM neutral massive 
gauge boson. Contributions that arise from Z2, Z§ and are subdominant or negligible. However, we keep 
their contributions in the formulation for completion. It also serves the purpose of showing that the addition of the 

6 In this analysis all three families of leptons are all localised with the same bulk mass parameter which also allows one the added 
benefit of escaping constraints from lepton FCNC. 
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higher KK modes of the neutral gauge bosons can only produce infinitesimal contributions and careful thought and 
taking into consideration the nature of the dependence of the contributions on the mass of the gauge states will 
yield that these contributions steadily decrease in magnitude as we move up the tower and hence fail to introduce 
considerable enhancements. The corrections from the higher fermion modes in the KK tower can also be similarly 
argued to be subleading. 

We feel that it is not only important to discuss the effects of the warped extras dimension in the rare decays 
of charm mesons but also correlate it to similar effects in other flavour sectors like beauty and strange. In the 
following sections we look at correlations between some observables that we derive in this work with observables 
in the AF = 1 processes in beauty and strange as derived in the works of Blanke et. al. [23]. particularly useful 
for this approach has been the generosity of the authors of these works in sharing their parameter sets which 
has allowed us to correlate observables in charm, beauty and strange at exact parameter points over the whole 
parameter space. The results are indeed quite interesting and will be discussed in the following sections. 

5.1 D° 77 

Although we do not have explicit calculations of the contributions of ND from the warped extra dimension to 
D° — > 77, we feel it is necessary for us to discuss what kind of effects can be expected in this channel for two 
reasons. 

• It is a process in which the models with a warped extra dimension can contribute only at the loop level and 
have no tree level contribution. 

• The SM LD contribution to D° —> {J, + /i~ is driven by the total branching fraction of D° —> 77 [49,50], 
and, as pointed out before in [52], is independent of whether the latter is generated by SM or by ND. Hence, 
enhancements to SM LD contributions to D° —> can be generated by ND contributions to the D° —> 77 
channel. 

Using the notations from [52], the contributions of this model can arise in the 1PR and the 2PR (or 1PI) 
contributions. The cited work clearly shows that the 1PR contribution dominates the 2PR contribution by about 
an order of magnitude and ND of the LHT-like framework can hardly enhance the 2PR over the 1PR. A similar 
argument holds here. Both the 1PR and the 2PR will get contributions only at the loop level and both from the 
same sources. 

• From the mixing of the heavy charged boson states W Y± and W ± with the W ± which gives the physical 
states. 

• From the heavier physical states W 2 and W 3 driving charged currents with SM fermions. 

The possibilities of large contribution of ND from the warped extra dimension in AC = 1 processes are solely 
because FCNC's are generated at the tree level itself. However, for this channel, due to the structure of the 
contributing diagrams, it is not possible for tree level ND to enhance the rate. The first contributions come only 
at the loop level as additional gauge states in the loop for the 2PR diagram and enhancements to the effective 7 
vertex for the 1PR diagram. 

Of course as stated earlier and shown in Table 1, the SM LD contribution dominates over the SM SD contribution 
by three orders of magnitude and if ND has to make its effect felt, through what is essentially SD contributions, 
it has to enhance the SD rates by more than three orders of magnitude over what the SM can produce and this is 
not possible for ND of the type being discussed in this paper. 

Hence, the verdict is clear. D° — > 77 remains dominated by SM LD contributions within the current framework 
and hence the SM LD rate for D° — > remains controlled by purely SM contributions to D° — > 77 even in the 

presence of ND from the warped extra dimension. 

5.2 D° y+yr 

We have seen that the SM SD contributions to this channel are extremely tiny of O(10 -19 ) and the SM LD 
contribution is almost O(10 6 ) larger. If ND has to make its effects felt in this channel through SD contributions 
then it has to enhance the SM SD rates by more that six orders of magnitude. Figure 3 clearly shows that this is 
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Figure 3: Histogram of the branching fraction of D° — > /i + /J~ . The lighter region lies below the estimate of the SM 
LD contribution. The counts in each bin represents the number of parameter points that contribute to the range 
of values of BR(D° — > /i + /i~) in that bin. 



possible in a large part of the parameter space. Enhancements of O(IC) 1 ) — O(10 2 ) are possible over SM LD rates 
which is denoted by the darker region of the histogram. 

The contributions to this channel are driven by the tree level V — A contributions to the four fermion vertex 
(uc)v~a(p + U>~)v—A which come primarily from the the mixing of the Z l and the Z' with the Z yielding the 
physical Z\ state. The contributions of the other neutral gauge bosons are negligible to up to a few percent. The 
new contributions to the V — A and V currents in the effective Hamiltonian can be parametrized through the 
following functions: 

, V D _ \ - 97{m)g-{m) 
av d _ \- 97{m)[g+{m)- g-{m)] 



where 

9sm 



n \ V 2 



\[%-K Sill 



2 , 



and the functions gf q parameterize the difference or sum between the left and right handed couplings to the neutral 
gauge bosons, ZqiL^qjL^ and Zl\ R l^ R , for qi = u,d and / = /i; this notation differs from [23], where differences 
in Al,_r are used to parametrize the g l functions. These functions are obtained from the effective unitary matrices 
defined in Equations (50) and (51) for quarks and Equations (45) and (47) for leptons. Thus, the partial decay of 
the D° meson gets the following contribution: 



G 



_> 



47r sin 



BR(£° -> = ^-^( -J—27T- ) flrn^mn^jl- 4-^Y°_ A Y°* A , (53) 



with Yy_ A including the Yq (x) , Y\ (x) and AYy_ A with the appropriate CKM factors for the first two. 

Current searches for this channel are done with the LHCb detector and branching fractions above the SM LD 
rates can be accessible in the near future. This channel can be searched in super flavour factories too, either in Belle 
II or in any of the proposed super tau-charm factories. As mentioned earlier, what is important for this channel 
is that D° — > 77 be measured with at least order of magnitude precision so that the SM LD contribution to this 
channel can be understood, something that is achievable in BESIII in the near future and is definitely within reach 
of the super flavour factories. Hence, this is a good channel to look for ND from the warped extra dimension even 
if no ND effects are seen in B s — > as we shall show later in Section 8. 
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The authors of the reference [77] claim that models with a warped extra dimension cannot leave any effects in 
D° —> because C^ has no role to play here. Actually, C±o, which is the only Wilson coefficient that matters 

in this channel, gets very large contribution from the new degrees of freedom, much larger than what CV (or Cg) 
can get and hence D° — s- is enhanced significantly. 

5.3 D -> X u uu 




-16 -14 -12 -10 -8 



Log 10 [BR(D -> X u vv)] 



Figure 4: Histogram of the branching fraction of D — > X u uv. 



The scenario for D — > X u vv is even more interesting. As we have mentioned before, the SM rate for this 
channel is driven by SD operators and a detailed analysis with light quark loops is not expected to bring about 
much enhancements. The new contributions to the V —A and V currents of this channel in the effective Hamiltonian 
can be parametrized through the following functions: 



AX- = 



E 

m 

E 



[9i (m)+gf(m)]g q (to) 
8g 2 SM M m 

\gf (to) + g+ (to)] [g+ (to) - g~ (to) 



(54) 



Therefore, in this channel, the partial width gets the following contribution: 



Br(D -> X u vv) 



G%m 5 D 



1927r 3 r D V47rsin 2 



Q 



7ir 



(55) 



where Xy_ A includes both the SM and ND contribution and Xy comes from ND only. 

It can be seen from Figure 4 that ND from the warped extra dimension can bring about enhancements of 
O(10 5 ) — O(10 6 ) over the SM rates in a large part of the parameter space. As for the case of the branching 
fraction of D° — > this is possible even when ND of such kind can only make negligible or no contributions 

to B s — > 

The source of such enhancement is again the mixing of the ND gauge states with the SM Z boson producing tree 
level FCNC. Like in the case of analogous kaons and beauty decays [23], here too there is vector FCNC generated 
in the ND which is absent in the SM. However, numerically, almost all the enhancement due to ND comes from 
the V - A FCNC. 

Such large enhancements as we see here can be accessible in super flavour factories. They might be out of reach 
for the LHCb as the branching fraction is relatively small and the final state comes with two neutrinos. Observation 
of this channel will provide a very clean and clear signature of ND. If signatures in this channel are found at the 
higher values of the branching fraction shown in Figure 4, a strong case can be made for tree level enhancement too, 
as the loop level enhancement seen in [54] is about one or two orders of magnitude larger than the SM estimates. 
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5.4 D XJ+l- 
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Figure 5: Histogram of the branching fraction of D X u l + l . 

The branching fraction of D — > X u l + l~ is not expected to have a significant enhancement, as the SM contribution 
is dominated by the photon penguin contribution where ND can hardly make its effect felt. The other contributions 
are negligible in comparison to the SM and need to be enhanced by several orders of magnitude in order to be 
comparable to the SM photon penguin contribution in Cg. As can be seen in Figure 5, the effects of ND in the 
branching fraction is quite small, 0(10%) in most of the parameter space, quite negligible as it is evident from 
the peak of the distribution lying on the SM value for the branching fraction. The additional contributions to Cg 
and Cio from tree-level neutral currents arising from the warped extra dimension are parametrized in the following 
way: 

'AY,? + AY D 

AC 9 = 

AC W = - (AY, 1 '- AY{[.. i ). 

where 



v r 2 ; v - A -H^z°+Az°_ A ) 

sin W 

rD a vD 



„ \qt (m) — q7 (m)] \qt (m) — q~(m)] 

AZ{? = - V iyi y 1 — yi - J liyq V ' — yq 1 (58) 

8g 2 SM M m 

While Cio is enhanced by orders of magnitude and Cg gets significant contributions too, these tree level contributions 
fail to compete with the photon penguin contribution in Cg coming from the SM. This is quite like what we saw in 
the LHT-like models [53] where loop level effects contributed to enhancing Cg, Cio and Cj{C%). In addition to this, 
the branching fraction in these channels are dominated by long distance effects by orders of magnitude. Hence, 
new dynamics have no chance of showing up in the branching fractions of these channels. 

The effects of models with a warped extra dimension in this channel was studied in [77] . They have missed out 
quite a few subtleties. Firstly, the largest effect that this kind of ND will have on this channel is not through its 
contribution to the dipole operators in c —¥ 1*7. In fact the dipole contributions to C7 and Cs hardly play a role in 
determining the size of the branching fraction of D — > X u l + l~ even in the SM. It is the photon penguin contribution 
in Cg which sets the stage for the SM SD contribution. As we have shown, the photon penguin contribution in Cg 
is not enhanced significantly by these models over SM values. The real enhancement is seen in Cio through tree 
level neutral currents which, however, can at most become comparable to the SM SD photon penguin contribution. 
Secondly, the branching fraction in this channel is not expected to be enhanced by ND. What is expected to be 
impacted by the presence of ND are the asymmetries as we show next. 
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6 Asymmetries in D — > X u l + l 




Figure 6: Histogram of the Asymmetries in D — > X u l + l . 

It is a difference story altogether for the impact of ND on asymmetries in the D — > X u l + l~ channel. We saw 
earlier that the SM signatures in the asymmetries are extremely tiny. From Figure 6, it is quite obvious that ND 
intervention in the asymmetries is sizable or even large - i.e., orders of magnitude more than what the SM can 
produce. 

The forward backward asymmetry A PB can be enhanced to even 0(5%) in some of the parameter space. This 
enhancement can be understood from the enhancement of C w as A FB depends directly on the magnitude of C w . 
Tree level contributions from the mixing of the new gauze states with the SM Z gauge boson bring about such 
enhancements and can be seen to be quantitatively more than the loop level enhancement that we had seen in the 
LHT-like models [53]. It should be noted that in the SM, and in the absence of CP violation (a good assumption 
within the SM for these channels), A PB for the conjugate channels should have opposite signs. However, in Figure 
6, it can be seen that the asymmetry for both the conjugate channels are of the same sign. This is a clear indication 
of the existence of CP asymmetry in these channels. 

The latter is made clear by the up to 0(1%) CP asymmetry, A^ P , that can be seen from the graph in Figure 
6. Although CP asymmetry remains small for quite an insignificant fraction of the parameter space, there are 
possibilities of measurable CP asymmetry in a significant fraction of the parameter space too. This is not surprising 
as there are new phases coming from the new mass mixing matrices of the fermions affecting the neutral currents 
driven by new gauge states. 

As expected from large A FB and sizeable A^ P , the A PB is enhanced by orders of magnitude in this scenario 
of ND. As seen in the LHT-like models [53], this asymmetry can be as high as O(10 4 %) and can be both positive 
and negative although it does show a tendency of being positive in a lot more of the parameter space. This can 
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be attributed to the prevalence of positive Ap B over negative ones in a large part of the parameter space, but can 
also arise from Ap B being larger in D than in D in most parts of the parameter space 7 . 

All three asymmetries provide unambiguous signatures of ND as SM values for these are very negligible. Probing 
these asymmetries should be on any menu that includes the experimental study of rare charm decays. To probe 
these asymmetries experiments have to be sensitive to at least the SD contribution to the branching fraction. While 
it is not necessary to discern the LD from the SD dynamics in the branching fraction itself unless ND changes the 
shape of the distribution of the differential branching fraction, this reach is necessary as the asymmetries are 
generated by short distance dynamics primarily. While LHCb might be sensitive to such small branching fractions, 
it might be complicated to extract the asymmetries from crowded background. The cleaner environment and full 
coverage of super flavour factories would be much better for such measurements. 

7 Correlations between different AC = 1 observables 




Log 10 [BR(D^X u / + r)] Log 10 |A° CP | Log 10 |/4 C FB I 



Figure 7: Correlations between D — s- X u l + l observables. The black dot in every graph represents the SM 
contributions to the corresponding observables. 

Having studied the observables themselves, it is instructive to take a look at the correlation between these 
different observables too, for three reasons: 

• Enhancement of one observable is not always accompanied by the the enhancement of another, even in the 
same decay mode. Hence correlations tell us a lot about what sets of observables to look for given a subset of 
the parameter space. 

• Correlations give us an insight into the underlying theoretical framework, in this case, the Wilson coefficients 
and their relative enhancements and dominance. Generalisation can follow from such an overview. 

• Most importantly, it is necessary to look into correlations between observables in different flavour sectors. 
This serves not only as a probe of flavour non-universality but also gives directions in search strategies in 
experiments. 

The first two we shall address in this section while the last case we shall study in the next section. 

7 This can be concluded only since a very tiny part of the parameter space has less than 1% 
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Amongst the decays studied in this work, the three body inclusive mode D — > X u l + l~ is the richest in contri- 
butions from different operators. Studying the correlations between the branching fraction and the asymmetries in 
this mode gives us an insight into how ND contributes to Wilson coefficients C-j{C%), Cg and Ciq. It is important to 
note here that the ND scenario under consideration in this work has new sources of FCNC at tree level only. Hence, 
significant contributions cannot be made to photon or gluon penguins unlike the case of LHT-like models [53] . The 
only significant contributions that are made are to Cg and Ciq. Cg is dominated by the SM photon penguin as 
mentioned earlier and hence the only Wilson coefficient that gets a significant enhancement over its SM value is 
Ciq. However, new phases can (and do) appear in both Cg and Cio due to the new fermion mixing matrices. 




Log 10 l^ c F Bl Log 10 |/\ C CP | Log 10 |A cp FB | 



Figure 8: Correlations between asymmetries in D —> X u l + l~ and BR(Z?° — > The region below the black 

line shaded in grey represent the part of the observable space where ND contribution to BR(_D° — > fi + n~) is 
subdominant to SM LD contribution to the same. 

It is clear from the tight correlation, in Figure 7, between the branching fraction in D — > X u l + l~ and the Ap B 
that the enhancements in both come from an enhancement of Cio as the latter is linearly dependent on C\g. As 
mentioned before, Cg is already large in SM and does not get enhanced much by ND. The correlations between 
Aq V and the other variables are very loose as the former is sensitive to the phases in Cg and Ciq. The only other 
set of observables that seem to be correlated are the Ap B and where one seems to be increasing with the other. 
The former is sensitive to the magnitude of C\g while the latter is sensitive to the phase in C\g which points at 
significant enhancements to both the magnitude and the phase in C\g from ND. 
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Log 10 [BR(D -» X u vv)] 



Figure 9: Correlations between BR(D — >• X u vv) and BR(D° — s- . The black dot represents the SM estimate 

of BK(D — > X u vv). The grey region represents the part of the parameter space where only BR(£> — > X u vv) gets 
enhanced above SM estimated while ND in BR(Z3° — > remains subdominant to SM. 

In Figure 8 we look at the correlation between the asymmetries in D — > X u l + l~ and D° — > /i + /J~. The 
grey band denotes the part of the parameter set that fails to overcome the SM LD contribution to D° — > fj, + fi~ . 
Therefore the points in the white region of each of the graphs represent enhancements to both the branching fraction 
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of D° —> and the asymmetries in D — > X u l + l~ that, if observed, would hint at the presence of ND. All the 

asymmetries show possibilities of enhancements to large values while the branching ratio of D° — > /i + fi~ remains 
large and distinguishable from SM LD contributions. The only two pairs of observables that are clearly correlated 
are the branching fraction of D° /i + fj,~ and Ap B since both depend purely on the size of Cio- 

The case for the the correlations between the branching fractions of D a —> /i + ' \i~ and D — > X u vv is shown in 
the plot on the right in Figure 9. In the grey region the branching fraction of D° — » lies shrouded in SM 

long distance dynamics. We see here that there are significant parts of the parameter space in which simultaneous 
enhancements to the branching fraction of both D° — > /i + /J~ and D — > X u vv are possible due to the tight correlation 
between the two observables. It can also be seen that even if ND shows up in D — > X u vv, it can fail to overcome 
the SM contribution to D° — » in a large part of the parameter space. However, it should be kept in mind 

that in the latter regime, it will be quite difficult for experiments to measure the branching fraction of D — > X u vv 
due to its small size. 



8 Correlations between strange, charm and beauty 

It has been long known that the intervention of ND in the different flavour sectors can be different. Firstly, there 
are possibilities of non-universal scenarios of flavour dynamics that can manifest themselves and leave different 
signatures in different flavour sectors. Secondly, manifestations of ND also need to dominate over SM contributions 
to flavour observables to be seen in flavour dynamics. While this is not possible in some flavour observable, some 
others, specially from charm dynamics, are ripe for such occurrences. A correlation study between AS" = 1 and 
AB = 1 observables was done in [23]. It is quite instructive to also compare flavour dynamics in two different 
sectors. 




Figure 10: Correlations between BR(B S — > ) and BK(D° — > ) . The black dot represents the SM LD 

estimate of BR(D° — > The y axis is well within the 95% C.L. values for the current LHCb results for 

BR{B S -> v+n-) [32]. 

In Figure 10 we show two observable and their correlations. The plot shows the correlation between BR(_B S — > 
and BR(£>° — > The values on the y axis lie well within the 95% C.L. values from the current 

experimental measurement of the branching fraction of B s — > announced recently by LHCb [32]. 

1.1 x 10~ 9 < BR(L>° ->■ /i+aO < 6.4 x 10~ 9 at 95% C.L. (59) 

There is not much correlation between the two observables. However, it is clear that even if the experimental 
errors are reduced and the branching fraction of B s —> is very close to the SM expectation and ND effects 

are indiscernible in that channel, large enhancements can show up in the branching fraction of D° — > and 
even dominate over the SM LD contributions to the same. This stems from the fact that SM contributions to 
B s — > (like in most of beauty dynamics) are quite large. Hence, ND can hardly make orders of magnitude 

enhancement there. The enhancements to B s — > ji + jjT can clearly be seen here to be of 0(10%) while in the 
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same parameter space it is orders of magnitude for charm. Hence, new dynamics can manifest itself in charm 
while keeping a low profile in beauty and vice versa, although the latter is a more difficult situation to deal with 
experimentally. We advocate the study of ND in charm lest nature choses the former. 




Figure 11: Correlations between BR(A'+ ->■ tt+vD), BR(A l ^vv) and BR(D° • The black dot 

represents the SM LD estimate of BR(D° — > /i + /i~) along with the SM estimates for BR(A + — > 7r + ^P) and 
BR(Kl —> tt q vv). The grey band on the left plot is the la experimental bounds on BR(A + — > ir + h>v). 

As mentioned earlier, a comparison between FCNC in strange and beauty dynamics was done in [23,24]. It 
was shown that the observables in the strange and beauty sectors are anti-correlated and can show enhancements 
in complementary parts of the parameter space. However, even in the AS — 1 observables studied in the koan 
sector, namely the branching fractions of Kl — > M + M~> Kl ~ * t^^vv and Kl — » ir°l + l~ the enhancements are 
modest and of at most an order of magnitude. The last two are compared in Figure 11 with BR(D° — > 
The experimental limits for these two koan decay modes are [89,90]: 

BR(A+ n+vv) = 1.7 ± 1.1 x 10" 10 (60) 
BR(A L ->■ n°vv) < 2.6 x 10~ 8 at 90% CL (61) 

to be compared to the SM values of [91] 

BR(A + — > ir + vv) = 7.83 ± 0.82 x 10~ n (62) 
BR(A L -> 7T°w) = 2.49 ± 0.39 x 10" 11 (63) 

While it is true that if the experimental value stays close to the central value for BR(A + — > tt+vD), this kind 
of ND cannot make large enhancements to BR(D° — > /i + /i _ ) it is also true in some parts of the parameter space, 
even if ND leaves very small signatures in the former, the latter can be enhanced by orders of magnitude in a large 
part of the parameter space. With experimental limits being quite close to the SM, the case of ND suffers the same 
fate as in beauty: there is a possibility of enhancements for ND but it lies shrouded in the shadows of the SM. 
We emphasise again, the case of charm is very different. Where ND can only leave modest effects in strange and 
beauty, it can leave a severe impact in charm, sometimes orders of magnitude beyond the reach of the SM. 

9 Analysis of the parameter space 

While we see enhancements in quite a few observables in rare charm decays, it might be instrumental to ask at 
this point whether such enhancements lead us to a preference for any particular part of the parameter space for 
the generation of observable new dynamics. In Section 4 we displayed the distribution of the Yukawa couplings 



24 



10 Inference 



c Q3o.3 



0.5 



0.4 



0.2 



0. 




0.50 



-0.55 



Figure 12: The bulk mass parameters of the three families of quarks and their correlation to enhancements in 
BR(D° — > n + fi~). The figure on the left shows the bulk mass parameters for the left handed quark doublets and 
the one on the right shows those for the right handed singlet quarks. The green (lighter) points represent the part 
of the parameter space that can generate enhancements to BR(D° — > beyond the SM LD estimates. 

and the bulk mass parameter that was used in this analysis. In Figure 12 we take another look at the space of the 
bulk mass parameters in the light of ND contributions to BR(_D° — > /i + /i~) beyond the SM LD estimates. 

The green (lighter) points represent the part of the parameter space in which we see enhancements to BR(D° — > 
/! + /i~) beyond the SM LD estimates. We clearly see, from the figure on the left, that there is a strong preference 
for the bulk mass parameter of the second family of left handed SM quark doublets to be at lower values close to 
0.5 when such enhancements occur, while there is no clear preferences for the values of the other two quark doublet 
bulk mass parameters. The figure on the right also shows that such enhancements prefers lower values of ct, the 
bulk mass parameter of the right handed top quark. 

The Yukawa couplings do not show any trends and the parameter points generating enhancements in BR(Z)° — > 
beyond the SM LD estimates are more or less randomly distributed in over the entire parameter space. 

10 Inference 

In this analysis we have looked at a particular flavour structure of possible ND from a model of warped extra 
dimension and tried to see what signatures it can leave in the dynamics of the charm quark. However, this analysis 
goes beyond just examining the latter effects in one particular model. In essence, this analysis can well achieve 
some model independence if we look into the details of the ND signatures left in the observables that we have 
studied. 

Let us first take a look at the flavour structure we have used. The model essentially has a set of new gauge 
bosons, not only from new gauge sectors but also from the appearance of KK towers. While it has both new neutral 
gauge bosons and charged gauge bosons, the former is more significant in an analysis of FCNC as such is generated 
at tree level itself, which is a novel signature of the ND absent in the SM. So it is sufficient to consider only the 
neutral gauge bosons in the analysis of FCNC. 

Although the fermion sector also has new degrees of freedom, none of that contribute to FCNC at any degree 
of significance as the tree level diagrams with SM fermions and new gauge bosons are dominant by far. While 
the actual size of the FCNC depends on multiple factors like the overlap of the fermionic profiles with the gauge 
boson profiles, the localisation of the Higgs, the Yukawa parameters and the KK mass scale, many of these are not 
independently scalable as there are multiple experimental constraints that tie these together. We shall go ahead 
with some generalisations of this flavour structure and its signatures in charm changing neutral currents (CCNC) 
based on two arguments, firstly, the nature of the spectrum of the degrees of freedom and secondly, the KK mass 
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scale. 

The localisation of the gauge bosons are fixed once the KK mass scale is fixed along with the boundary conditions 
to match up with the SM gauge content and, in the case of the current model, a wish to preserve custodial symmetry 
in the IR brane. We shall not consider the localisation of the Higgs as small changes in the same within allowed 
model constraints do not have any large effects in the observables studied 8 . The localisation of the fermions can 
change within certain limits as can the Yukawa parameters. The last two parameter sets are what determine the 
size of CCNC, the KK mass scale remaining the same. 

From the above arguments, we can define a generic flavour structure: 

• New neutral and charged gauge boson states. 

• Overlap of the SM fermions with the new gauge boson states generating tree level FCNC. 

• New angles and phases coming from fermion mixing matrices producing flavour non-diagonal neutral currents. 9 

We claim that as long as models with a warped extra dimension are built of this flavour structure it will contribute 
to CCNC significantly. Inclusion of more neutral gauge fields can only enhance the current effects seen but not by 
orders of magnitude as long as such inclusion follows all the requisite constraints from electroweak precision tests 
and other precisely measured flavour observable. Here we have considered a non-ad-hoc flavour structure that is 
not "tuned" to give effects in the up type quarks. Moreover, we have used a parameter set whose effects have 
already been studied in [23,24]. Hence it gives us a comprehensive picture of the effects of these kinds of flavour 
structures in the dynamics of both up type and down type quarks. 

Our second argument brings us to consider the limitations of the LHC in direct searches. One motivation for 
keeping ND at the TeV scale is to make sure that we see it at the LHC direct searches. In many models with a 
warped extra dimension it is possible to keep low KK scales using different theoretical technologies. Yet, one is led 
to wonder what will happen if we fail to discover these low KK scales, i.e., nature has "heavier" plans for us. 

The most simplistic way one can argue the KK scale dependence in the observables we have studied is to state 
that the amplitudes depend roughly on the inverse square of the KK mass scale and hence the branching fractions 
and the asymmetries depend on the same to the fourth power. The KK mass scale that we have used is about 2.45 
TeV. It would be quite justifiable to say that if the KK scale really lies at around 5 TeV or more 10 , direct searches 
at LHC will start having a problem seeing new degrees of freedom. Yet the effects in charm dynamics would at 
best be lowered by an order of magnitude. Considering what the numbers tell us, charm dynamics would still be 
in the game for showing ND effects even if direct searches at LHC and possibly, both beauty and kaon dynamics 
would be out of the game within the ambits of such a ND scenario. Suffice it to say that charm has been both 
cursed and blessed by tiny SM signatures. 

11 Summary 

The possibilities of the existence of ND coming from a warped extra dimension is an intriguing solution to many 
problems and questions that face us today. Yet these models are severely constrained by what has already been 
measured. While efforts are being made to both see the new degrees of freedom in direct searches and probe new 
dynamics in flavour dynamics, we put in our share of trying to see the impacts of these structures on the dynamics 
of the charm quark and their features. 

We also, in general, strongly disagree with the conclusion in [77] that the "only" place the models with a 
warped extra dimension can show its effect are in CP asymmetries in D — > X u j besides AAcp 11 . We think we 
have convincingly argued in this paper that models of this kind can have large effects well beyond the SM estimates 
in multiple rare decays including D° — » D — » X u vv and asymmetries in D — > X u l + l~ which can be within 

experimental reach in the near future. Also, we propose a study of higher multiplicity hadronic decays of the 
charmed meson both experimentally and theoretically, although a lot of tools need to be developed for these. 

In this work we have shown that: 

8 For a study of the dependence of CP violation in D° — > tt + tt~ / K + K~ on the Higgs localisation cf. [77]. The fact that the Higgs 
localisation does not matter much in the effects that we see in our work shows that these effects are independent of what needs to be 
done to generate AA^p . 

9 In the absence of any new angles and phases, the model reduces to the MFV kind which we do not study. 
10 We are assuming the reality of the existence of a warped extra dimension. 
11 For detailed studies of the connection between these two observables cf. [92,93]. 
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• These kind of flavour structure can leave large effects in charm, sometimes orders of magnitude larger than 
what the SM can generate. 

• While these models leave moderate effects in beauty and strange dynamics, the effects in charm need not be 
moderate. This can be achieved even without giving the up type quark sector a special dynamical advantage. 

• The effects in charm dynamics are not tightly correlated with those in beauty and strange, i.e., even when 
this kind of ND can leave negligible contributions to both beauty and strange dynamics, it can leave very 
large contributions to charm dynamics. 

• These tree level effects coming from this class of models can be larger than the loop-level enhancement that 
we saw in the LHT-like models for the same observables. 

An important part of this work is the study of correlations between the dynamics of the down type quark sector, 
i.e., B and K decays and the dynamics of an up type quark in D decays. Not only do we connect dynamics of 
different quarks through such an analysis, but we also connect the two sectors of quarks within a parameter space 
of new dynamics that obeys experimental constraints both from EWPD and flavour observables. We show that 
even while ND might not manifest itself prominently in the dynamics of the down type quarks, it can well bring 
about orders of magnitude enhancements to dynamics of the charm meson as we have highlighted above. It should 
be noted that the class of model we have studied here is not motivated to produce effects in any of the flavour 
sectors, collectively or singly. In spite of this, we not only see possible large impacts in flavour dynamics, but also 
a clear distinction between effects in different flavour sectors. 

While it continues to be true that accessing rare charm dynamics is statistically challenging and theoretically not 
well understood even at the level of the SM contributions, we have made some strong cases for ND intervention way 
beyond what the SM can generate. These effects not only lead to the observables coming within reach of current 
experiments and future super flavour factories, but also sidestep the problem of determining SM contributions 
theoretically due to their large size. 

Finally we want to state that indirect evidences for ND are based on flavour dynamics beyond the SM in 
general. However we do not like to go to a shopping mall to find anything that is just beyond the SM flavour 
dynamics; we greatly prefer to think about flavour dynamics that originates from a motivation to find solutions to 
the hierarchy problem or challenge of the SM. Previously we have worked about LHT to deal with models of one 
class, now we have done the same for the Warped Extra Dimension, a la, the class of Randall- Sundrum models. 
Now we shall wait and hope that nature conspires on our side. 
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